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ABSTRACT 

The quantity of dust in a spiral disk can be estimated using the dust's typical 
emission or the extinction of a known source. In this paper, we compare two 
techniques, one based on emission and one on absorption, applied on sections of 
fourteen disk galaxies. The two measurements reflect, respectively the average 
and apparent optical depth of a disk section. Hence, they depend differently on 
the average number and optical depth of ISM structures in the disk. 

The small scale geometry of the cold ISM is critical for accurate models of 
the overall energy budget of spiral disks. ISM geometry, relative contributions of 
different stellar populations and dust emissivity are all free parameters in galaxy 
Spectral Energy Distribution (SED) models; they are also sometimes degenerate, 
depending on wavelength coverage. Our aim is to constrain typical ISM geometry. 

The apparent optical depth measurement comes from the number of distant 
galaxies seen in HST images through the foreground disk, calibrated with the 
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"Synthetic Field Method" (SFM). We discuss what can be learned from the 
SFM measurement alone regarding ISM geometry. 

We measure the IR flux in images from the Spitzer Infrared Nearby Galaxy 
Survey in the same section of the disk that was covered by HST. A physical model 
of the dust is fit to the SED to estimate the dust surface density, mean temper- 
ature, and brightness in these disk sections. The surface density is subsequently 
converted into the average optical depth estimate. 

The two measurements generally agree and the SED in order model finds 
a mostly cold dust (T < 25 K.). The ratios between the measured average and 
apparent optical depths of the disk sections imply optically thin (r c = 0.4) clouds 
in these disks. Optically thick disks, are likely to have more than a single cloud 
along the line-of-sight. 

Subject headings: (ISM:) dust, extinction, ISM: structure, galaxies: ISM, galax- 
ies: spiral, infrared: galaxies, infrared: ISM 



1. Introduction 



The dust content of a spiral galaxy disk can be mapped either by the characteristic 
dust emission in the far-infrared (FIR) and sub-mm regimes, or by using the attenuation of 
known background sources. Both techniques have seen recent significant improvements in 
accuracy and sensitivity, with complementary results shedding light on the dusty interstellar 
medium in spiral disks. 

The emission from the interstellar dust in the disks of spiral galaxies has been charac- 
terized with increasing accuracy by several infrared space missions (IRAS, ISO and, recently, 
Spitzer), as well as the sub-mm observations of SCUBA on the JCMT. The improvements 
in spatial resolution and wavelength coverage have led to significant insight into the tem- 
perature components of the dust in spiral disks, and into the relation between dusty clouds 
and star-formation. The FIR emission from spiral galaxies has revealed that the dust can 
be described by two dominant thermal components: warm (25 K < T < 100 K) and cold (T 
< 25 K). Both the warm and cold components can be found in spiral arms, and there is a 
smooth disk of cold dust b etween these arms. Most of the dust mass in the spiral disk is cold 
(see t he review articles by lGenzel &: Cesarskyil2000i : iTuffs &: Popescull2005l : iPopescu fc Tuffs 
2005h . 



FIR and sub -mm observations of galaxies find indications (e.g., iTrewhella et al.l |2000| ; 
Alton et al.l Il998l ) or direct evidence of cold dust disks extending beyond the stellar disk 
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(iNelson et all 1199a iDavies et al.lll999l: IPopescu fc Tuffsl 120031) . The studies of edge-on spi- 
rals by iRadovich et al.l (120011 ) and IXilouris et al.l (119991 ) quantified the radial profile as a 
scalelength of the dust that is 40% larger than that of the starlight. The contribution of the 
cold ISM (T < 25 K) to the overall emission of spiral disks has been difficult to constrain be- 
cause of the degeneracy between dust temperature and mass. Hence, the cold ISM's relation 
with HI and their relative distribution remain unknown. 

Spitzer observatio ns, mostly from the Spitzer Infrared Nearby Galaxy Survey (SINGS, 



Kennicutt et al.l 120031 ). have already contributed greatly to the understanding of spiral 



disks. The relations between the tracers of cold dust (70 and 160 micron emission) and 
star-formation, both obscured (24 micron) and unobscured (UV and Ha emission) have 
already been studied with this multi- wavelength survey in several canonical g alaxies and 



their substructure: the starburst M51 ( Calzetti et al. 20051; Thorn 



design spiral M81 (I Gordon et al. 



2004 



Perez- Gonzalez et al 



faegan et al .ll2004h and M31 JGordon et alJbood) . the superwind in M82 fenge 



20061), and the dwarf NGC 55 (lEngelbracht et al. 



ey et alJbood ). the grand 



2006|), the rings of NGC 7331 



20041 ). bale et all (120051. 120071) discuss the 



bracht et al. 



SED of all SINGS galaxies over all available wavelengths. iDraine et al.l (120071 ) find ample ev- 
idence for dust in all the SINGS galaxies, with the gas-to-dust ratio related to the metallicity. 
They find no evidence for very cold (T < 10K) dust, however. 

These studies find ample evidence of cold dust throughout the optical disks of spi- 
rals but, interestingly, also outside them in various places: on the e dge of the optical 



disk (IThornley et al.ll2006l : iGordon et al.ll2006), outside M82's superwind ( lEngelbracht et al. 
20061 ) and extending beyond the stellar disk (IHinz et al.ll2006l ). 



Parallel to these investigations of dust emission has gone an observational effort to quan 
tify th e absorption by dust in spiral disks using known background sources. IWhite fc Keel 
(119921 ) proposed using occulting galaxy pairs for this purpose. Nearby occulting galaxy pairs 



were initially investigated with ground-based data, both images (jAndredakis & yan der Kruit 



1992 



Berlind et al. 1997 ; Domingue et al.fl999 ; White et aliboOO ) and spectra ( Domingue et al 



20001 ). Subsequently, with the Hubbl e Space Telescope (H S T), a more detailed pic ture of 
dust in these nearby disks emerged ( Keel Sz White 2001al lbl; Elmegreen et al. 2001 ). The 
results of these studies are that extinction is gray 1] when measured over disk sections greater 
than 100 pc, but resembles the Galactic extinction law at smaller scales - those that can 
only be resolved with HST. Arms are found to be more opaque than the general disk, and 



1 Gray extinction is equal attenuation at all wavelengths: there is no relation between color and measured 
optical depth. A color measurement is dominated by the lines-of-sight with the least extinction, while the 
independent extinction measure is dominated by those with the most extinction. In cases where many 
lines-of-sight are mixed, gray extinction is mimicked. 
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some evidence suggests that the dust disk is a fractal, similar to the HI disk. 



Gonzalez et al.l (119981 ) investigated the use of the calibrated number of distant galaxies 
seen though the foreground disk in HST images. The calibrate d counts of distant galaxies 



have been explored fur ther in the pr evious papers in this series (IGonzalez et al.l 1 19981 . 12003 



Holwerda et al. 2005a. b 



20061 ). Both the occulting galaxy t echnique and counts o f 



distant objects yield very similar opacities for disks and spiral arms (IHolwerda et al.ll2005bl ). 



In recent years, models have been developed to explain the Spectral Energy Distri- 



bution (SED) of ec 


goon spirals spanning wave 


lengths ranging from the UV to the FIR. 


(e.g., 


Popescu et al. 


2000; 


Misiriotis et al. 2001: 


Popescu et al. 200( 


); 


Vlisiriotis et al. 


2001; 


Tuffs et al. 2004: Boissier et al. 2004: 


Dasvra e 


t al. 2005: r 


jiffs ct 


al 


. 12004: Dasvra et al. 


2005|; 


Calzetti et al. 


2005; 


Dopita et al. 


2006b a; 


Draine & Li 


2007: 


Draine et al. 2007 


)• 



Three scenarios have been proposed to explain the discrepancy between the apparent ab- 
sorption in UV and optical wavelengths and the emission of dust in the FIR and sub-mm 
regimes: 



1. 



A young stellar po p ulatio n is embedded in the de nse plane of the disk . This is proposed 
by iPopescu et al.l (120001 ) and corroborated by iDriver et al.l (120071 ). The embedded 
young stars pump the FIR emission radiated by the dust plane. 

A strongly clumped dusty medium. The clumpin g would lead to underestimate the dust 
mass from optical extinction in e dge-on systems (IBianchi et al.ll2000bl ; IWitt fc Gordon 



UV to FIR SED (IBianchi et al 



2000l ; iMisiriotis fc Bianchil 12002); the dust mass would also be underestimated by a 



2000af ) 



3. A different emissivity of the cold dust grains -higher than canonical- in the FIR and 



sub-mm. A change of emi ssivity has been proposed for denser ISM regions (lAlton et al. 



2004; 


Dasvra et al. 


(Bendo et al. 


2006) 



20051 1 or, alternatively, for the lower density regions of the disk 



In all three of these scenarios, the dumpiness of the dusty ISM is an important factor. 

While in so me models the large scale structure of the dusty ISM has been somewha t 
constrained (e.g.. lXilouris et al.lll999l ; ISeth et al.ll2005l ; lBianchill2007l ; iKamphuis et al.ll2007l ). 
the small-scale geometry ("dumpiness") of the cold ISM remains unknown. Therefore, an 
estimate of the prevalent dusty cloud size for spiral disks would provide a constraint for the 
SED models of spiral disks. 

Given that SED and extinction techniques are sensitive to different characteristics of 
the dusty clouds in the spiral disk, a comparison between the optical depth derived from 
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these two methods has the potential to figure out the structure of the dusty ISM. Here, 
we compare the J-band optical depths for a section of the spiral disk, one derived from a 
SED model of the Spitzer fluxes ("average"), and one determined from the number of distant 
galaxies found in an HST image ("apparent"). The term "average" refers to what the optical 
depth proportional to the dust mass, uniformly distributed over the disk section. The term 
"apparent" means the effective optica l depth of the disk sect ion for a background uniform 



light source (original definitions from iNatta fc Panagial Il984l ) . The term "opacity" is used 



throughout our previous papers for the apparent optical depth measured over a section of 
the disk for its whole height. 

This paper is structured as follows: Sj2] discusses the Spitzer and HST data used. In §31 
the two different methods to derive optical depths are presented. We discuss the relation 
between dust geometry and galaxy counts in §U In $5j we present the derived optical depths; 
£0 presents a simple geometric model to interpret the results, and $7] lists our conclusions 
and future work. 



Data 



The dat a for this paper co me from two archives, the HST archive and the fourth 



data release flS 



Kennicutt et al 



NGS team! 120061 ) of the Spitzer Infrared Nearby Galaxy Survey (SINGS 



2003f . The r e is an overlap of fourteen galaxies between the SINGS sample 



and that of iHolwerda et all (12005fi). Two of these, N GC 3621 and NGC 5194, have two 
WFPC2 expos ures analyzed injHolwerda et all (ho05bh . The HST/WFPC2 data reduction 
is described in IHolwerda et al.l (j2005al ). The reasoning behind the HST sample selection 
from the archive is explained in §3.11 

The overlap between HST and Spitzer data is illustrated in Figured! with the WFPC2 
footprint projected on the 24 micron Spitzer images. Only the solid angle covered by the 
WF chips is used for further analysis (the PC chip is excluded). 

The Infrared Array Camera (IRAC) mosaic is made with the custom SINGSdither script, 
by M. Regan, that combines the s can images into a single mosaic using the "drizzle" algo- 
rithm (Fourth Data Release Notes. ISINGS teamll2006l ). The Multiband Imaging Photometer 
for Spitzer (MIPS) dat a products are calibra t ed, sky-subt r acted mosaics in all t hree b ands, 
reduced as described in iGordon et al.l (120051 ) ; iBendo et al.l (120061 ) ; ISINGS team! (120061 ) . The 
basic instrument parameters, pixel scale, and adopted PSF FWHM for the seven main Spitzer 



2 http:/ /sings. stsci.edu 
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imaging modes are summarized in Table HJ 



3. Analysis 

Two parallel estimates of the optical depth of disks are used in this paper: first, the 
apparent optical depth of the spiral disks is determined from the number of distant galaxies 
identified in the HST/WFPC2 images, calibrated with the "Synthetic Field Method" (SFM). 
Secondly, the optical depth of the same section of the spiral disks is derived from the dust 
surfa ce density, w h ich is a result of the SED model fit to the Spitzer fluxes using the model 



from lDraine fc Lj (120071 ). 



3.1. Galaxy counts: "Synthetic Field Method" 



In principle, the number of distant galaxies seen through a spiral disk is a function of 
the dust extinction, as well as the crowding and confusion by the foreground disk. Initial 
applicat ions of the number of distant ga l axies as an extinction tracer were on the Magellanic 
Clouds Jshaplevl Il95ll : IvVesselinkl Il96ll : Irlodge! ll974J : iMacGillivravl Il975h. but thev lacked 
accuracy. The "Synthetic Field Method" (SFM) was developed by I Gonzalez et al~ ( 1998 ) 
to correct an extinction measurement based on the number of distant galaxies in an HST 
image, for the effects of crowding and confusion by the foreground spiral disk. 

The "Synthetic Field Method" follows a series of steps. First, the number of distant 
galaxies in an HST science field is determined. Second, a series of simulated ("synthetic") 
fields are made. In each of these fields, a typical background (e.g., the Hubble Deep Field) 
is first dimmed by a gray screen and added to the science field. Third, the added distant 
objects are identified in these synthetic fields. The fourth step is to measure the relationship 
between the number of these identified synthetic distant galaxies and background dimming. 
From this relation and the original number of actual distant galaxies found in the science 
field, an average opacity can be found. It is important to remake the synthetic fields for each 
science field because the crowding and confusion issues are unique in each case. 

An additional uncertainty in the resulting average extinction measurement is the cos- 
mic variance in the intrinsic number of distant galaxies behind the foreground disk. The 
uncertainty due to cosmic variance can be estimated from the two-point correlation function 
of distant galaxies and folded into the Poissonian error. The cosmic variance uncertainty 
is of the same order as the Poisson statistica l error for small numb ers (For a complete dis- 
cussion on the uncertainties of the SFM, see lHolwerda et al.ll2005al ). Therefore, single-field 
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SFM measurements remain uncertain, but a meaningful conclusion can be drawn from a 
combined set of science fields. 



We have applied this method successfully on archival WFPC2 data. iHolwerda et al. 



( 12005bl ) present the average radial opacity profile of spiral disks and the effect of spiral arms. 
The spiral arms are more opaque a nd show a strong radia l dependence, while the more 
transparent disk shows a flat profile. IHolwerda et al.l (12005a ) compare HI radial profiles to 
the opacity ones, and conclude that no good relation between disk o pacity and HI surface 
density radi al profiles can be found . However, IHolwerda et al.l (12005a ) find that the sub-mm 
profile fromlMeiierink et al.l (120051 ) generally agrees with their opacity measurements of M51. 
Holwerda et al.l (j2005el ) compare the relation between surface brightness and disk opacity; 
this relation is strong in the spiral arms, but weak in the rest of the disk. 



Gonzalez et al.l (120031 ) predicted, based on simulated data, that the optimum distance 



for the application of the SFM with current HST instruments is approximately that of 
Virgo. The identification of background galaxies suffers in closer disks, as the stellar disk 
becomes more resolved, compounding confusion. This optimum distance, combined with 
the availability of deep HST/WFPC2 images from the Cepheid Ex t ragalac ti c Distance Scale 



Key Pr oject, resulted in the samp le presented in IHolwerda et al.l (l2005bl ). IHolwerda et al. 
(2005d) confirmed the results from Gonzalez et al. ((2003 ) using this sample, with foreground 



disks spanning distances between 3.5 and 35 Mpc. A selection effect of the Key Project is 
that the majority of the HST science fields is concentrated on spiral arms and exclude the 
centers of the galaxies. 

Here, we present average extinction values for the whole WFPC2 field-of-view, mini- 
mizing the uncertainties to the extent possible. In our initial papers, we did not apply an 
inclination correction to the optical de pths because the corre ction depends strongly on the 
dust geometry (see the discussion in IHolwerda et al.l l2005bl ). However, in this paper we 
assume a simple dust model in §61 The appropriate inclination correction -x cos(i)- has 
been applied to the points in Figures H] and [51 and in Table HI The uncertainties in the tables 
and figures reflect the 1-sigma confidence levels produced by the combination of the Poisson 
error and the cosmic variance of background galaxies. 
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3.2. SED optical depth estimate 



The average disk optical depth is derived from the Spitzer observations. First, the 
surface brightnesses within the WF chips' footprint are measured^] (see Figured]). Second, 
these are converted into a dust surface density using an SED model. Third, this surface 
density is translated into an J-band optical depth. 

All the IRAC and MIPS data are convolved to the poorest resolution of the 160 micron 
observations (see Table [TJ, and the pixel scale is set to 9". This is done with the gauss, 
wcsmap and geotran tasks, under IRAF. Subsequently, the overall flux is measured in the 
WFPC2 field of view (Figured]). Because the L-shaped aperture is a highly unusual one, the 
aperture correction remains uncertain, but not negligible, since the FWHM at 160 micron 
(40'.'0) is of the order of the aperture diamet er (3 x 1.3 x 1'.3 in an L-shape). Published 
aperture corrections for the IRAC instruments (IHora et al.ll2004j ) overestimate the correction 
for e xtended object s (jjarrett 20051 ) . Here we use the aperture cor r ection s for extended sources 
from ljarrettl (hoosh for IRAC fluxesfl and from iMuzerolle et al.l (hoosh for the MIPS fluxefl 
(see also Table [[]). Table [2] gives the average surface brightnesses for the seven Spitzer 
channels in the field-of-view of the three WF chips of the WFPC2 array. The uncertainties 
are derived from the variance in the s ky. Generally, the surface brightnesses agree with the 
results presented by iDale et al.l (120051 ) for the entire disks. 



The second step is to convert these surface brightnesses to a dust surface density. Ini- 
tially, we fi tted only the MIPS fluxes with two blackbodies and derived surface densities 
from these (IHolwerda et al.ll2006l ). However, a more rigoro us treatment of the IR fluxes can 
be done with a SED model, such as the one presented in iLi &: Draind ( 120011 ) . This model 
uses the physics of grain heating and reradiation, and a model distribution of grain sizes 
and types. The updated version from lDraine fc Li] (120071 ) has been fit to the data, and the 
results are shown in Figure [2] The derived stellar and dust surface brightnesses, dust surface 
density, and mean temperature are presented in Table [3] Dust surface densities are between 
0.1 and 1.4 x 1O 6 M kpc~ 2 , with mean temperatures between 14.6 and 17.8 K. 

The mean dust te mperatures are obtained from the mean radiation scaling, U, in 
the iDraine fc Li] (120071 ) model (T2 in their equation 18). The model uses a distribution 



3 The PC part of the WFPC2 array is not used in either analysis. 
4 http:/ /spider. ipac.caltech.edu/staff/jarrett/irac/calibration/index. html 



3 http :/ /ssc . spitzer . caltech.edu / mips /apercorr / 



6 The lDraine fc Lil (]2007h model uses a distribution of scaling values (U) of the local Interstellar Radiation 
Field to calculate the irradiation that the grains see. U is the average of this scaling distribution. 
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of temperatures and grain sizes. Therefore the mean temperature is an indication of the 
thermal equilibrium point of the bulk of the dust. Most of the dust is cold (T < 25 K). 

These results are an obvious improvement over a simple single-temperature fit, but the 
model parameters in Table [3] are still not fully constrained; disagreement between data and 
model at the PAH peak at 8 micron could be an effect of metallicity or the presence of a 
bright HII region. The 70 and 160 /i m fluxes hint at colder or more dust in the disk (Figure 
[2]). There are three caveats to the fits: (1) a lack of sub-mm data, (2) the single color ISRF 
used and (3) averaging over different types of emission regions, i.e., HII regions and the 
general disk. 

Additional sub-mm data of comparable quality, needed to better constrain the model 
and especially the cold dust emission, will not be available until SCUBA2 starts operations 
on the JCMT and the launch of Hersche§ At presen t, comparable-qu ality sub-mm map s 
are available only f or NGC 5194 jMeiierink et al.ll2005h and NGC 7331 faegan et al.ll2004h . 
Draine et al.l (120071 ) discuss SED models with and without sub-mm data. 



The second caveat in the derivation of dust mass is the assumption of a constant color 
for the interstellar radiation field (ISRF) illuminating the emitting dust. In reality, dust 
grains deeper in a dust structure will encounter radiation field that is not only dimmed but 
also reddened and, hence, will contribute less flux to the FIR emission. Locally, the ISRF 
will also depend on the age of the nearby stellar population. Due to the reddening of the 
ISRF deeper in the cloud, dense clouds could contain more dust mass in their centers than 
inferred from just the FIR emission. A dditional sub - mm o bservations will help resolve this 
uncertainty in the SED optical depth. iDraine et al.l (120071 ) discuss fits of the model to the 
total fluxes of the SINGS galaxies with and without additional sub-mm data. They find that 
the FIR estimate underestimates of 1.5 times the dust mass more (5 out of 17 cases, notably 
NGC 3627 and 7331 of our sample) than overestimates it (only M51 of the 17). 

The third caveat is that the model values are an average over many different types of 
ISM regions, each with a different heating mechanism, dust structure and comp osition (e.g., 
photo -dissociation regions, cirrus, and star- forming regions in spiral arms). The lDraine &: Li 
( 120071 ) model's assumptions hold better for some regions than for others, but we use the 
results as "typical" for these disks. The relative contribution of PAH emission to the SEP 



is a f unction of ISM ge ometry as well as irradiation and composition (e.g., ISilva et al. 



19981 : IPiovan et al1l2006h . Together with a better-constrained FIR/sub-mm S ED one could 



constrain ISM geometry solely from the relative contributions to the SED. IDraine et al. 



7 Herschel data will be especially valuable will be especially valuable as it will not suffer from night sky 
structure, which is of similar angular size as these disks. 



-10- 



( 120071 ) discuss the application of the lDraine fc Lil (120071 ) model to whole disks of the SINGS 
galaxies. 

The dust surface density is translated into an average opt ical dep t h wit h the absorption 
cross section per unit dust mass, k Q 6 S (A), and grain albedo from lDraind (120031 ) for the Johnson 
J-band (865.5 A.): r m = K a b s x (1 — Albedo) x Mdustl 'area. These optical depth values are 
presented in Table HI 



4. Cloud size and the SFM 

It continues to be difficult to constrain dusty cloud geometry from models of either 
extinction or emission. In this section, we review what can be learned, solely from the SFM 
measurements, about the geometry of the extincting medium. 



In iHolwerda et al.1 (l2005bl ). two indications that the dust disk is clumpy are identified: 

(1) the average color of the distant galaxies is independent of the disk opacity implied by 
their number, and (2) the measurement of disk opacity is independent of inclination. The 
lack of a relation between the average opacity of the disk and the average color of the 
detected galaxies can be explained by two scenarios: (1) some of the background galaxies 
are completely blocked by large clouds, and some are not. The color measurement is done on 
background objects that do not suffer from extinction and, hence, reddening. Alternatively, 

(2) all background galaxies are dimmed by clouds smaller than the projected distant galaxies. 
Consequently, some of the distant galaxies are dimmed enough to drop below the detection 
threshold. Any detected galaxy's color is, however, measured from mostly unreddened fluxjfl 
However, we note that a relation between the reddening and derived extinction from the 
distant galaxies is difficult to detect because (1) the spread in colors of distant galaxies is 
substantial and (2) color is measured from the detected galaxies -automatically the least 
dimmed- whereas opacity is measured form the missing galaxies. 

The inclin ation effect on the appa rent optical depth from the number of distant galaxies 



is minimal (see IHolwerda et al.ll2005bl . , §5.1 and Figure 3). Assuming a thin layer of optically 
thick clouds, the apparent optical depth of the disk, measured from the number of distant 
galaxies, is dominated by the apparent filling factor of clouds. The projected filling factor 



8 The effect of cloud geometry on the luminosity function of the detected distant galaxies is, in principle, 
another possible way to distinguish between cloud geometries, but this is complicated by the fact that the 
detection limit of distant galaxies is dominated by the field p roperties (brightness an d granularity), rather 
than by the d imming by the dust. See also the discussions in IHolwerda et all (|2005dh . and Appendix B in 
Holwerdal (fcoosh . 
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does not change much with inclination: a flat cloud covering 40% of a certain disk section still 
covers 40% of the inclined section. Only when the height of the cloud becomes important 
-when the inclination is closer to edge-on-, does the apparent filling factor change. This 
explanation for the lack of an inclination effect in the opacity profiles does not depend on 
the size of the clouds. It could be a single, large, cloud or many sm all ones in the plane of the 
disk. However, the optical depth values in lHolwerda et al.l (j2005bl ) are from different sections 
of the disks -although generally centered on a spiral arm- and the effect of small inclination 
differences coul d well h ave been masked by comparing different regions in the disks. In 
Holwerda et al.l (l2005bl ) we did not apply an inclination correction because it depends on 
the assumed dust geometry. In this paper we do assume a dust geometry and hence make 
an inclination correction (§6]). 

The simulations in the SFM assume a gray screen, an uniform unclumped dust layer 
with opacity equal in the V a nd I bands . The SFM opacity measurements in this paper are 
based on such simulations. In Holwerda ( 2005 ). we ran a series of simulations on NGC 1365^1 
to characterize the effect of average cloud cross-section on the number of distant galaxies 
observable through a disk. 

Figure [3] shows the effect of cloud size, expressed in pixels, on the simulated relation 
between average opacity and number of distant objects. In each simulation, we fix a cloud 
size and vary their number to increase disk opacity. An ensemble of unresolved clouds is 
effectively the gray screen. For clouds resolved with HST, more than 2 pixels^], the relation 
between opacity (cloud filling factor) and number of distant objects becomes much shallower. 
The same number of distant galaxies observed would then imply a much higher opacity 
of the disk. Because the SFM (calibrat ed with a gray screen) generally agrees well with 
measurements from overlapping galaxies (IHolwerda et al.ll2005bl ). it seems unlikely that the 
disk's opacity is predominantly due to large, resolved clouds. A pixel of 0'.'05 at the distance 
of NGC 1365 (18 Mpc) is 4 pc in linear size. It is therefore implied that the structure of 
the ISM responsible for the disk opacity measured with the SFM varies in optical depth on 
scales of ~10 pc or less. 

From the SFM measurements alone, the cloud geometry is impossible to determine. 
Only when additional information is used -e.g., the general agreement with the occulting 
galaxy technique- it favors small (unresolved) scales for the clouds. Therefore, to constrain 
cloud geometry, information from two different techniques needs to be combined. 



9 NGC 1365 is one of the more distant galaxies in the Holwerda et al. ( 2005hh sample, and we ran these 
simulations as a validation of the gray screen synthetic fields. 

10 Pixel scale is 0'.'05 in our data, after the drizzle reduction. 
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5. Optical depths 



Table H] presents the optical depth estima tes from the Synthetic Field Method and the 
SED model JPraine! liool braine fcU l2007h . for the WFPC2 field-of-view. The optical 
depths range between 0.1 and 3.5 magnitudes in the /-band. The measurements are for 
different parts of different spiral disks (Figure [TJ, explaining in part the range in values. 

The optical depth estimates presented here may appear high for the Johnson /-band, 
compared with other extinction estimates (e.g., those from inclination effects or reddening), 
but these are (1) for the entire height of the disk, and (2) generally centered on a spiral arm. 
Typical extinction values in the /- band are several te nths of a magnitude for a dust screen in 
front of the stellar spiral disk (e.g.. lMeyer et al.ll2006l ). Two of the derived SFM opacities are 



negative, possibly the effect of an overdensity of distant galaxies behind the target galaxy. 
The discrepancy in NGC5194-1 may be due to misidentification of background galaxies, as 
they are difficult to identify in this field. 

Figures |4] and [5] show the values of disk opacity by both methods, over the same section 
of the disk. Both methods generally agree within the uncertainties of the measurements. 
T he agreement is b etter than our initial estimate from a blackbody fit to the MIPS fluxes 



m 



Holwerda et al.l (120061 ). The general agreement and the mean temperature of the dust 



(Table [3]) imply that most of a disk's opacity is due to the cold dust in the disk. 



6. Model of cloud geometry 

The relation between the apparent and average /-band optical depth measurements -the 
first from the number of distant galaxies in HST images and the second derived from the 
Spitzer SED- could reveal the nature of the prevalent structure in the ISM. 

The a priori assumptions are that (1) all dust structure is transparent to the FIR 
emission from which the dust surface density is estimated in the SED m odel, and (2) the 



entire volume of the cloud emits in the FIR. We adopt Model C from iNatta &: Panagia 



( 119841 ). in which a randomly distributed series of clumps covers the area. These authors 
define two optical depths: (1) the typical optical depth (r m ), i.e., the average optical depth 
that is proportional to the dust mass; and (2) the apparent optical depth (f), or the optical 
depth if a uniform layer would cover the area. 

Our two measurements of optical depth -SED and SFM- correspond to these two optical 
depths, average and apparent. An optical depth based on the SED depends on the dust mass 
within the area, and hence corresponds to r m . The SFM optical depth is the apparent optical 
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depth averaged over the area, and hence r. These two optical depths need not be the same, 
and their relation is an indication of how clumped the medium is. 

Let us assume a number of small dust structures with a height (h), an average grain 
cross-section (a), and a grain emissivity (Q). The grain number density in the clouds is 
denoted by rid, and the average number of clouds in a line-of-sight is n. We assume all 
clumps have the same optical depth r c : 

t c = h n d a Q. (1) 
The average optical depth (r m ) is then: 

r m = n x r c , (2) 

and the apparent optical depth can be derived if one assumes a Gaussian distribution of the 
number of clouds along the po ssible lines-of-sight and sum the contributions of all clouds 



(see equations 15 and 17-19 in iNatta fc Panagialll984l ): 



f = n x (l-e~ Tc ), (3) 
with asymptotic values for both optically thin and thick clouds: 

r c — > 0, => r m -» 0, f — >• T m , (4) 

r c -> oo, => r m -> oo, f n. (5) 

We note that the inferred dust mass in the disk {M^ust) is proportional to the number of 
clouds (n) and the cloud optical depth (r c ). The ratio of the apparent over the average 
optical depth is: 

f 1 — e~ Tc 

f = ( 6 ) 

leaving only the optical depth of the clouds (r c ), and hence the grain density {n d ) and the 
cloud size (h) as the variables. Our fit to the relation between SFM and SED optical depth 
estimates only has r c as the variable (Figure H]). 
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We can now use equation [6] to derive r c from a fit to the optical depths from SED (r m ), 
and from the number of distant galaxies (f). We want to answer three questions. Are the 
clouds in the disks typically optically thin or thick? If all disks are equal, what is the implied 
cloud optical depth? How many clouds lie typically along a given line-of-sight? 



6.1. Optically thick or thin clouds? 

Optically thin clouds (r c << 1) result in a ratio of optical depths close to unity: 

f l-e- Tc 1 - (1 - T r 



i; (7) 



optically thick clouds (r c >> 1) result in a ratio of: 



r l - e ,c 
— = < I. 



Figures H] and [5] show how a majority of the data exhibit a ratio of order unity. The cold 
dust mass in the SED model could be better constrained with additional sub-mm information 



(jDraine et al.ll2007l ). However, the ratios are, within the errors, consistent with optically thin 



clouds in most of the disks. 



6.2. Cloud size 

Figure H] shows the fit to the ratios of apparent to average optical depths, with r c as 
the single fit parameter, as per equation [61 For simplicity, we assume here that all disks are 
made up of similar clouds, and that there is no difference between arm and disk regions. The 
negative SFM measurements are excluded from the fit. The best fit is for a cloud optical 
depth of r c =0.4 with, on average, 2.6 clouds along the line-of-sight (if the two negative 
points are included, the values change, respectively, to r c = 0.56 and n = 1.9). The inferred 
value of t c is likely to be a mean between higher values in spiral arms and much lower values 
in the disk. Optically thick disks have more clouds along the line of sight, while optically 
thin disks harbor a single cloud. 

Reasonable values for the parameters in equation [1] are: iij k 5 x 10~ 3 grains m~ 3 , 
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a = 0.03 /mi 2 and Q = ^ (±§ f = 1.5 x 10~ 3 with (3 = 2 jHildebrandlll983l lR. The value 
of 0.4 for t c implies a cloud height, h, of ~ 60 pc! Much larger clouds could be resolved 
in extinction maps of these disks based on stellar reddening. In the case of NGC 3627, 
NGC 5194, NGC 6946 and NGC 7331, there is a clear spiral arm in the reddening map 
(IMeyer et al.l 120061 ; iHolwerda et al.l 120071 ) ; the other reddening maps are smooth or do not 
extend out to cover the whole WFPC2 pointing. 

A more likely scenario is that there is an inverse relation between the cloud density rid 
and s cale h. Such a relatio n can be seen in giant molecular clouds (GMC) of our own Galaxy 



(e.g, iSolomon et al.lll987l ). In this case, the single optically thin cloud can be replaced by 
smaller optically thick ones. We note that the largest GMCs are of the order of 60 pc. 

The value of 60 pc clouds appears in contradiction to the implied size in the SFM cali- 
bration (~10 pc). However, the typical cloud size can be 60 pc and still the disk opacity can 
change over smaller scales, if several partially overlapping clouds are seen in projection. The 
inverse relation between scale and density would also help make the two scales compatible. 

Our data are consistent with optically thin clouds (r c = 0.4) and their average opacity 
value implies a typical cloud size that is unresolved with Spitzer in our galaxies. 



6.3. Cloud Numbers 

The above fit to the relation between f and r m indicates that, on average, more than 
one cloud is needed along the line-of-sight in most disks (n = 2.6). However, the assumption 
was that r c had a single value for all the disks. In §6.11 we argued that the ratio between 
f and r m implied the r c is optically thin. It logically follows that optically thick disks must 
have more than a single cloud along the line-of-sight. 

Figure [5] illustrates the effect of number of clouds along the line-of-sight (n) on the 
relation between f and r m , when r c is freely increased from to r m /n (the lines have been 
drawn according to equations [2] and [3]) . We note that all disks are consistent with many 
clouds along the line-of-sight -including the optically thin ones, also because they lie in the 
optically thin cloud (f/r m ~ 1, r c < 1) regime. 



11 There is su bstantial disc u ssion in the literature about the value of f3 (e.g., 
fixed at 2 in the lDraine fc Lil (|2007l ) models. 



Bendo et al 



2003h . P was 
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6.4. Potential Improvements 

There are many refinements to be made to the simple model presented here. Some 
improvements for future comparisons between these two measurements of optical depth are: 
(1) a distribution of cloud sizes for both the SFM calibration, as well as in the model 
explaining the ratio between SFM and SED optical depths. A model d istribution can be 



taken from obser vations of GMCs in our own Galaxy and nearby ones (IHeyer et al.1 12001 



Rosolowskyi 120051 ) . The cross-section distribution could be used in SFM measurements in 
the future. A cross-section distribution can only be applied, if the foreground disk is at a 
single, fixed distance; only one counts though a single foreground galaxy are used. There are 
three face-on spirals with enough solid angle in HST imaging as well as additional Spitzer 
data: M51, M81 and M101. (2) The SED model can be much better constrained with 
additional sub-mm observations, the opportunities for which will expand dramatically in the 
near future (SCUBA2 on the JCMT and Herschel satellite). (3) The effects of grand design 
spiral arms and galactic radius could be identified in a single disk; the comparison SFM 
and SED is not made for different sections of the disks combined. (4) A future SED model 
can take into account the reddening of the interstellar radiation field, as it penetrates the 
ISM. This would require a comprehensive treatment of the ISM structure in addition to its 
temperature, composition and irradiation. 



7. Conclusions 

To constrain models of the spiral disk's energy budget with typical values for the size 
of dusty clouds in the ISM, we compare two techniques to extract the average and apparent 
optical depths of a section of spiral disk. From the comparison between SFM and SED 
results, we conclude the following: 

1. The SFM's calibration alone implies projected cloud scales predominantly unresolved 
with HST (of the order of 10 pc in NGC 1365, see §Uand Figure [3]). 

2. The dust responsible for the disk's opacity is predominantly cold (T < 25 K, Table EJ). 

3. The average and apparent optical depths of these disk sections, measured from SED 
and SFM respectively, generally agree ( Figure H] and ED- This implies generally optically 
thin clouds (r c < 1, AO) . 

4. The fit to the ratio between apparent and average optical depth measurements, f/r m , 
indicates a cloud optical depth, r c , of 0.4, more than a single cloud along the line-of- 
sight, and a cloud size of ~ 60 pc. If several partially overlapping clouds are seen in 
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projection through the disk, the disk's opacity will change over smaller scales, consistent 
with conclusion 1. 

6. Optically thick disks appear to have more than a single cloud along the line-of-sight 
(Figure [5]) and optically thin disks may have several clouds as well. 



Future work using counts of distant galaxies through a foreground disk could be used to 
find cold dust structures at larger galactic radii, provided a sufficiently large solid angle has 
been imag ed with HST/ACS's superb resolution^ Notably, the ACS data on M51, M81 and 
M101 are very promising for such an analysiso Spitzer observations of these nearby disks 
are also available, making a similar comparison between SED and apparent optical depth 
possible for portions of these disks. The typical cloud scale for spiral arms or disk sections or 
as a function of galactic radius could then be found. The SCUBA-2 instrument has recently 
been installed on the James Clerk Maxwell Telescope. A project with SCUBA-2 to map 
the SINGS galaxies in two sub-mm bands will improve future SED modelling of these spiral 
disks significantly over the SED models presented here. 
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Fig. 1. — The footprints of the WPFC2 camera on board HST, on the 24 micron images from 
the MIPS detector onboard Spitzer. Most of the HST images do not include the center, and 
are pointed on a spiral arm. NGC 3621 and NGC 5195 (M51) have two separate WFPC2 
fields associated with them. The PC chip, i.e, the small chip in the nook of the "L" of the 
three WF chips, is not used for the SFM analysis, nor as part of our Spitzer aperture. 
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Fig. 2. — The Spectral Energy Distribution in the 7 Spitzer bands (I RAC/MIPS) , for ea ch 
of the WFPC2 apertures in Figure [TJ The best fitting model from iDraine fc Lil (120071 ) is 
shown. The relevant parameters for each fit are in Table [31 
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Fig. 3. — Different synthetic fields with different dust disk models from iHolwerdal f 20051 ) . A 
set of synthetic fields is made using dimmed HDF images. This dimming can be a smooth 
uniform gray screen (triangles), or some distribution of opaque clouds (curves). The relation 
between the number of distant galaxies from the HDF that can still be retrieved and the 
average opacity of the dimming depends on the assumed model. Large, resolved clouds block 
fewer background objects given the same filling factor. As a result, one a higher optical depth 
is implied by the intersection of the curve and number of observed background galaxies. The 
distance to NGC1365 is 18 Mpc, so a pixel of 0'.'05 is equivalent to 4 pc in linear size. The 
scales in the above simulations correspond, therefore, to clouds with cross-sections with a 
radius of 2, 4, 8, 100 and 400 pc, respectively 
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Fig. 4. — The mean and apparent optical depth in the /-band, from SED (r m ) and SFM 
(f). The best fit with equation [6] to the ratios of these values is also shown (r c = 0.4). The 
average value for number of clouds along the line-of-sight, n, is then 2.6. If the negative 
SFM values are included in the fit, the cloud optical depth rises to r c = 0.6 and the average 
n becomes 1.9. 
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Fig. 5. — Model values of r m and f, from equations [2] and El for fixed values of the average 
number of clouds along the line-of-sight (n = 0.5, 1, 2, 3, 4). The cloud optical depth r c is 
left to vary proportionally to r m (r m = n x r c , equation [2]). Hence maximum cloud optical 
depths are 6, 3, 1.5, 1 and 0.75 respectively. If optically thin clouds are made up of optically 
thin clouds, it seems necessary that more than one cloud will lie along the line-of-sight. 
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Table 1. SINGS data description. IRAC and MIPS pixel scales, PSF FWHM, and 
aperture corrections (the factor by which the fluxes are multiplied). 



Instrument 


pixel 


PSF 


aperture 


(band) 


scale 


(FWHM) 


correction 




i 


2 


3 


IRAC (3.6) 


0"75 


2"5 


0.9 


IRAC (4.5) 


0"75 


2"5 


0.9 


IRAC (5.8) 


0"75 


2"5 


0.7 


IRAC (8.0) 


0"75 


2"5 


0.75 


MIPS (24) 


1"5 


ero 


1.16 


MIPS (70) 


4"5 


1870 


1.2 


MIPS (160) 


9"0 


40'/0 


1.4 



1 Pixel scales were set by the SINGS team. 

2 FWHM values for the IRAC are conserva- 
tive estimates. Actual FWHM values are better 
than 2'/5. 

3 The IRAC values are the initial results from 
T. Jarrett (private communication), but they 
do not differ substantially from the final results. 
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Table 2. Spitzer channel surface brightness in the WFPC2 aperture (3 WFs xl.3' x 1.3'). 



Galaxy 


3.6 /im 


4.5 /im 


5.8 fim 


8 (im 


24 fim 


70 fim 


160 fim 




MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 


MJy/sr 



NGC 


925 





17 


± 





01 





19 


± 





08 





12 


± 





05 


3 


79 


± 


2 


65 





33 


± 





06 


4.76 ± 0.44 


21.92 ± 0.91 


NGC 


2841 





39 


± 





04 





33 


± 





09 





88 


± 





59 


3 


42 


± 


2 


13 





16 


± 





04 


4.54 ± 0.34 


34.27 ± 0.58 


NGC 


3031 





56 


± 





03 





39 


± 





06 





25 


± 





07 


2 


07 


± 


1 


11 





23 


± 





04 


3.74 ± 0.37 


23.73 ± 1.38 


NGC 


3198 





13 


± 





02 





12 


± 





04 





11 


± 





04 


2 


28 


± 


1 


20 





32 


± 





05 


3.48 ± 0.36 


20.09 ± 0.46 


NGC 


3351 





55 


± 





05 





15 


± 





16 





33 


± 





07 


5 


81 


± 


4 


11 


1 


11 


± 





06 


12.54 ± 0.54 


43.55 ± 0.68 


NGC 


3621-1 





24 


± 





01 





12 


± 





01 





40 


± 





05 


1 


27 


± 





37 





70 


± 





05 


9.68 ± 0.46 


41.38 ± 0.87 


NGC 


3621-2 





15 


± 





01 





17 


± 





01 





24 


± 





05 


1 


08 


± 





37 





66 


± 





05 


11.04 ± 0.46 


36.81 ± 0.87 


NGC 


3627 





97 


± 





04 





68 


± 





06 


1 


20 


± 





06 


6 


72 


± 


2 


89 


5 


80 


± 





06 


49.58 ± 0.61 


152.63 ± 0.95 


NGC 


4321 





55 


± 





02 





39 


± 





04 





64 


± 





06 


4 


45 


± 


2 


23 


2 


62 


± 





05 


25.12 ± 0.46 


100.53 ± 0.75 


NGC 


4536 





11 


± 





04 





15 


± 





07 





72 


± 


3 


45 


3 


64 


± 


2 


89 





19 


± 





05 


2.35 ± 0.52 


10.97 ± 0.60 


NGC 


4559 





05 


± 





03 





01 


± 





01 





06 


± 





05 


1 


51 


± 


1 


05 





10 


± 





05 


2.65 ± 0.38 


11.80 ± 0.53 


NGC 


4725 





11 


± 





02 





09 


± 





05 





28 


± 





20 


1 


76 


± 


1 


29 





16 


± 





05 


1.59 ± 0.35 


12.50 ± 0.49 


NGC 


5194-1 





52 


± 





01 





32 


± 





01 





80 


± 





05 


2 


61 


± 





22 


2 


64 


± 





04 


28.61 ± 0.38 


124.68 ± 1.36 


NGC 


5194-2 





56 


± 





01 





34 


± 





01 





88 


± 





05 


2 


76 


± 





22 


3 


26 


± 





04 


36.35 ± 0.38 


144.30 ± 1.36 


NGC 


6946 





43 


± 





03 





27 


± 





01 





71 


± 





09 


2 


31 


± 





48 


2 


92 


± 





07 


28.88 ± 0.53 


80.83 ± 2.13 


NGC 


7331 





14 


± 





01 





10 


± 





01 





10 


± 





05 


1 


39 


± 





88 





23 


± 





04 


3.49 ± 0.67 


17.05 ± 0.79 



Table 3. Model output: st ellar and dust surfa ce brightness, dust surface density, dust 
mean temperature (T 2 in lDraine fc Lil 120071 ). fit quality, the parameters of stellar 
irradiation -minimum and maximum of the distribution, mean irradiative field, fraction of 
dust exposed to more than U m i n . Models can be found at 
[http: / / www.astro.princeton.edu/~draine/ dust /irem. html 



Galaxy L s tar/area L dust /area M dust /area T \ 2 U m in Umax U 7 model 

name (Lq/^pc 2 ) (Lq/^pc 2 ) (M^/hpc 2 ) (K) name 

xlO 8 xlO 8 xlO 6 xlO 6 xl0~ 3 



NGC 0925 


0.8 


0.2 


0.2 


15.7 


2.8 





70 


10 


0.8 


5.8 


U0.70.1e7.MW3.1.20 


NGC 2841 


1.8 


0.3 


0.5 


14.6 


4.15 





70 


10 


0.8 


1.0 


U0.70.1e7.MW3.1.20 


NGC 3031 


2.7 


0.2 


0.3 


14.6 


2.32 





50 


1.0 


0.5 


1.5 


U0.50.1e6.MW3.1.30 


NGC 3198 


0.6 


0.2 


0.3 


15.0 


3.73 





50 


10 


0.6 


7.5 


U0.50.1e7.MW3.1.30 


NGC 3351 


2.7 


0.5 


0.4 


16.2 


1.9 





80 


1.0 


1.0 


1.6 


U0.80_le6.MW3.1.10 


NGC 3621-1 


0.7 


0.5 


0.4 


16.2 


2.45 


1 


00 


1.0 


1.0 


2.6 


U1.00.1c6.MW3.1.60 


NGC 3621-2 


0.7 


0.4 


0.2 


16.9 


10.91 


1 


20 


10 


1.3 


7.3 


U1.20.1c7.MW3.1.30 


NGC 3627 


3.7 


2.1 


1.1 


17.0 


2.11 


1 


00 


10 


1.4 


24.5 


U1.00.1e7.MW3.1.30 


NGC 4321 


2.2 


1.1 


1.0 


16.0 


1.91 





70 


10 


0.9 


15.6 


U0.70.1e7.MW3.1.30 


NGC 4536 


0.5 


0.1 


0.1 


15.7 


2.32 





80 


10 


0.8 


2.7 


U0.80.1e7.MW3.1.60 


NGC 4559 


0.2 


0.1 


0.1 


15.7 


1.76 





80 


1.0 


0.8 


0.0 


U0.80.1e6.MW3.1.30 


NGC 4725 


0.5 


0.1 


0.2 


14.6 


3.01 





50 


1.0 


0.5 


0.0 


U0.50.1e6.MW3.1.60 


NGC 5194-1 


1.6 


1.4 


1.2 


16.0 


2.33 





80 


1.0 


0.9 


9.9 


U0.80_le6.MW3.1.50 


NGC 5194-2 


1.8 


1.7 


1.4 


16.0 


1.83 





80 


1.0 


0.9 


11.7 


U0.80.1e6.MW3.1.40 


NGC 6946 


1.4 


1.2 


0.5 


17.8 


0.65 


1 


50 


1.0 


1.9 


20.5 


U1.50.1e6.MW3.1.50 


NGC 7331 


0.6 


0.2 


0.2 


15.4 


2.15 





70 


10 


0.7 


4.5 


U0.70.1e7.MW3.1.30 
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Table 4. The apparent (t) and average (r m ) optical depths in /-band measured in the 
WFPC2 field (3 WFs 1'3 x 1(3), uncorrected and corrected for inclination. 



Galaxy 


T 

(SFM) 


T 

xcos(i) 


7~m 

(SED) 


7"m 

xcos(i 


T/r m 

) 


NGC 0925 


-0 4 3 


-0.2 


0.6 


0.3 


-0.63 


NGC 2841 


0.7°4. 4 


0.4 


1.4 


0.7 


0.53 


NGC 3031 


o.8°4. 6 


0.4 


0.9 


0.5 


0.88 


NGC 3198 


7 3 

u - ' -0.3 


0.3 


0.7 


0.4 


0.97 


NGC 3351 


i - 1 -0.5 


1.1 


1.0 


1.0 


1.06 


NGC 3621-1 


z - u -0.6 


1.1 


1.0 


0.5 


2.05 


NGC 3621-2 


± " u — 0.3 


0.5 


0.7 


0.4 


1.38 


NGC 3627 


i qo.6 


1.0 


3.1 


1.7 


0.61 


NGC 4321 


2 2°- 7 


1.7 


2.8 


2.2 


0.76 


NGC 4536 


o.8°4. 3 


0.4 


0.3 


0.1 


2.58 


NGC 4559 


o.i°4. 3 


0.1 


0.3 


0.1 


0.39 


NGC 4725 


n 7 - 3 

u -'-0.3 


0.5 


0.5 


0.4 


1.38 


NGC 5194-1 


-0 4 - 4 


-0.4 


3.3 


3.0 


-0.12 


NGC 5194-1 


1 3 05 


1.2 


3.9 


3.5 


0.34 


NGC 6946 


1 o - 5 

i - u -0.5 


0.8 


1.3 


1.0 


0.76 


NGC 7331 


u -' 3 -0.3 


0.2 


0.5 


0.3 


0.62 



